Magnetization structure of a Bloch point singularity 
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Switching of magnetic vortex cores involves a topological transition characterized by the presence 
of a magnetization singularity, a point where the magnetization vanishes (Bloch point). We ana- 
lytically derive the shape of the Bloch point that is an extremum of the free energy with exchange, 
dipole and the Landau terms for the determination of the local value of the magnetization modulus. 
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The interest in the dynamics of magnetic vortices was 
renewed by the discovery of fast core reversal by a vary- 
ing external excitation (magnetic fielcP or spin currentpl) . 
Micromagnetic simulations of vortex core switching^ re- 
vealed that the underlying mechanism, the annihilation 
of a vortex-antivortex pair^, needs the mediation of a 
magnetization singularity: a magnetic monopole or Bloch 
pointPThis magnetization structure was first studied by 
Feldtkellerpl who showed that it is mainly determined 
by the exchange energy; later on Doring^ calculated this 
specific energy and showed that its value would be a 
topologically invariant. He considered a family of mag- 
netization textures differing in their local rotation angle 
7 (with respect to the radial direction) and found that 
minimization of the demagnetization energy density se- 
lected a specific angle 7 « 112°. However, any approach 
within the micromagnetic approximation (the magneti- 
zation strength is at its saturation value) cannot account 
for the internal structure of the singularity that imposes 
the vanishing of the magnetization vector. In order to 
investigate the region near the singular point Galkina et 
alP included the Landau magnetic energy, although they 
neglected the demagnetization term, and showed that the 
magnetization vector modulus increases linearly with the 
radial distance from the origin. Therefore, to under- 
stand the topological transitions between different vor- 
tex states, for which magnetic monopoles are required, 9 
it is important to go beyond the micromagnetic approx- 
imation. In this paper we compute the magnetization 
field of a Bloch point taking into account the exchange, 
Landau and demagnetizing energies. We obtain two so- 
lutions, the first one, corresponding to a local minimum 
of the energy density, is characterized by a linear mag- 
netization modulus near the origin and by an essentially 
azimuthal magnetization configuration with a rotation 
angle 7 (incidentally rather close to the one found by 
Doring-^); the second one, also linear near the center, 
but with a radial magnetization vector (hedgehog Bloch 
point), is valid over a finite spherical region, and corre- 
sponds to a local maximum of the energy density. 



and of the magnetic potential <f>, 
'A 



dV 



W + / L (M) 



Fh, (1) 



where dV is the volume element, A is the exchange energy 
constant, /l is the Landau energy density 
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a = a(T) is in general a function of the temperature, a < 
in the relevant ferromagnetic state, 6 is a dimensional 
constant, and the energy of the demagnetizing field 



F H = T H [M,®} 



-Mo 



dV (M ■ H + — 



(3) 



with H(r) = — V$(r) the magnetic field. We introduce 
the following units: length, £ = (A/hq) 1 ' 2 ; magnetiza- 
tion, M s = (-a/26) 1 / 2 ; and energy, E = (-aA/2b)£. In 
this units system the free energy becomes, 
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where f(M) = -M 2 + M 4 /2, and v = |a|//z is the 
only nondimensional parameter of the system; it can be 
written as 

v = £ 2 /£l £ = {A/^) 1 ' 2 , £ = {-A/a) 1 ' 2 , (5) 

where £ is related to the exchange length £a (£a = v^£) 
and £q is the characteristic length of the magnetization 
intensity M = \M\ variation, as will be demonstrated 
below. 

The equilibrium distributions of the magnetic poten- 
tial and the magnetization field are determined by the 
variation of T with respect to and M . The varia- 
tional derivative of Q with respect to leads to the 
Maxwell equations, 



In order to determine the magnetization field M(r) of 
a Bloch point in a ferromagnetic nanostructure we con- 
sider the free energy T — J-[M, $] as a functional of M 



V 2 $ = V • M, (6) 
in the magnetic domain, and, at the surface boundary 
n ■ M — n ■ AH, (7) 
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where h is the normal and AH the discontinuity of the 
magnetic field. The variational derivative of Q with 
respect to M leads to 
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This equation can be transformed into an integro- 
differential equation for the magnetization field using the 
explicit solution of (6][7) in terms of M, 
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where prime variables refer to the magnetic domain V 
and its boundary dV . 

The problem now is to determine the structure of the 
Bloch point as a particular solution of Eq. This 
can be done by introducing an appropriate ansatz. To 
compute the demagnetizing field we consider a simple 
geometry, we take for V a sphere of radius R; in this re- 
gion, we choose a magnetization field that generalizes the 
Feldtkellei^ and Doring^l ansatz, adding a r-dependent 
magnetization modulus, 



M = M(r)m 7 (0,0) 



(10) 



where the unit vector m 7 has the topology of a Bloch 
point and satisfies the condition to be an extremum of the 
exchange free energy (neglecting other terms in J 7 ). The 
simplest one-parameter solution, depending on a rotation 
angle 7, can be written asP 

m 1 — (cos(0 + 7) sin 6*, sin(0 + 7) sin#, cos 9), 

in cartesian coordinates (x,y,z), where 9, <fi are, simul- 
taneously, spherical angles for both r and m 7 . In the 
following calculation we use the notations 

c±(7) = -(l±coS7), c 2 (7) = ^(l + 2cos 7 ) 

and represent vectors in spherical coordinates (f,6,4>): 

c + (7) + c - (7) cos 26* 



-c_(7) sin 29 
sin 7 sin 9 



(11) 



where 7 = corresponds to the hedgehog configuration, 
and 7 = arccos(— 11/29) to the spiral Bloch point found 
by DoringP 

The magnetization intensity vanishes near the singu- 
larity and reaches its saturation value far from the origin 
(at a distance r ^> £q/£). Therefore, a solution of |8]) 
must satisfy: 



lim M{r) -> 0, lim M(r) M s (i/), 
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where the saturation magnetization M s is in general a 
function of v. Inserting the ansatz (10 1 into d8l one ob- 



tains the integro-differential equation, 
DM + 2v (M - M 3 ) 



H = 0, 
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where 
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is the radial part of the laplacian with angular momen- 
tum 1 = 1. 

To find the demagnetizing potential we insert the 



ansatz (10 1 into the equation for the potential (|9j) , and 
note that the magnetization volume V • M and surface 
r ■ M charges, can be expressed as an expansion over the 
spherical harmonics Y§(9 : (f)) and Y2(9,<fi): 

$(r) = 2Vnc 2 ( 7 )F(r)Y ° + ^ c _( 7 )G(r)r 2 ° (15) 

where we define the functions of the radial coordinate r, 
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(We use the standard notation 
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= r if 

r < r' , etc.) From this formula one immediately deduces 
the demagnetizing field, 

c 2 (j)F'(r) + c_( 7 )(i + cos29)G'(r) 
H = I _? c _( 7 )G(r) sin 29 | (IS) 



(in a spherical frame), that can be used in Eq. (13 1 to 
obtain the following four equations: 
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Different magnetization textures can be solution of these 
equations depending on their characteristic length scales 
and possessing different energies. The existence of the 
parameter v (usually large for ferromagnetic materials) 
allows the separation of three regions: a singular core 
region r < io/(, characterized by a rapid variation of 
the magnetization modulus; an intermediate region (the 



3 



0.5 



I 


L 








J 


2 4 6 

(2v) 1/2 r 


8 10 




FIG. 1. (Color online) Bloch point for the minimun energy 
structure, 7 / 0. (top) Magnetization intensity, in units of 
Ms as a function of the stretched coordinate \f2vr in units of 
£, and its angular distribution (inset) for 7 = 113°. (bottom) 
Nondimensional demagnetization energy as a function of the 
rotation angle 7. 



micromagnetic core) £q/£ < r < 1; and an external re- 
gion r > 1, dominated by the dipolar energy. In the 
singular core region (r < 1), where one can formally take 
v — > 00 , the demagnetizing terms are negligible compared 
to the exchange and Landau ones. It is then possible to 
distinguish between two cases: (i) sin 7 7^ 0, implying 
0(M) — 0, that leads to a local solution, valid near the 
center of the Bloch point, and, as we shall see, corre- 
sponding to a minimum of the energy; and (ii) sin 7 = 0, 
that allows for a global solution, corresponding to a max- 
imum of the energy. 

First we consider case (i), the minimum energy local 
solution. From (19 1, the assumption sin 7 7^ 0, implies 



O(M) — 0, and compatibility with the other equations is 
possible in the singular region, where the demagnetizing 
field is small. We have thus, to find a solution of 

0{M) = DM(r) + 2v[M(r) - M{rf] = 0, (24) 



with the boundary conditions 112} . A simple scaling 
transformation r — > [2v) 1 / 2 r allows to scale out the 
parameter v. This implies that the function M — 
M((2v) l l 2 r) is universal, which can be considered as a 
generalization of the invariance statement by Doring 7 . 
Let us consider first the singular region r < £q/£ <C 1 
where one can assume that M(r) = C(y)r is a linear 
function of the radial coordinate. This choice is moti- 
vated by the property that a linear magnetization satis- 
fies DCr — 0, and is then an asymptotic solution of (24 1. 



The constant C depends trivially on the parameter v\ us- 
ing the above radial scaling, one finds Oy) = {2v) 1 / 2 C 00 
where the universal constant must be determined as 
an eigenvalue of 0(M) = satisfying (12). It is im- 



portant to note that a linear magnetization amplitude 
near r = implies G ~ F ~ r 2 , meaning a posteri- 



ori that the demagnetizing terms in ( 20] 22 1 are indeed 
negligible compared to the exchange and Landau terms 
in the singular region. We consider now the core region 
r < 1, where, in accordance with (12 1, the modulus M 



approaches the saturation value M = 1, canceling the 
Landau term in (24 1. We remark that for the core re- 
gion r w 1, the large v condition ensures that M reaches 
its saturation value and that the demagnetization terms 
are also negligible. Indeed, the characteristic length for 
which the magnetization saturates may be estimated by 
(2^)!/ 2 r « O(l), or r « O(£ /£) that usually is very 
small. Therefore, the solution of (24 1 is consistent with 
the system (19J22) throughout the core region. 

The universal magnetization profile, a solution of 
0(M) — satisfying (12 1 with M s — > 1, was computed 
numerically and is shown in Fig. [T] We observe that the 
magnetization saturates for values of {2v) x / 2 r w 4; as a 
consequence the validity condition of the solution is well 
verified for v > 10. The numerical profile is also consis- 
tent with the linear magnetization near the center with 
Coo = 0.506, and a saturation M = 1 for large r. It is 
worth noticing that in this case the value of 7 is not de- 
termined, showing that the ansatz (JTo|) represents a one- 
parameter family of solutions for the internal structure 
of the Bloch point. However, the demagnetizing field, 
although negligible for the determination of the magne- 
tization radial profile, should select a specific value of 7 
in order to minimize the free energyP The relevant part 
of the free energy is the demagnetizing field density (the 
Landau and exchange terms are independent of 7), 

— (M • H + H 2 /2)dV = (H 2 /2)dV. (25) 

In the inner region, where the magnetization is given by 

M(r) = {2 V ) 1 / 2 C 00 r, (26) 

the demagnetizing field writes 



H = Cr 



-ca(7)+c-(7)(!-§oos20) 
|c_(7) sin 29 




(27) 



After integration over the spherical angles, the part of 
the energy density depending explicitly on 7, is 



6(7) = (28cos7 + 18cos2 7 + 29)/75, 



(28) 



(see Fig. [T] bottom panel) whose minimization gives 7 = 
arccos(— 7/18) s» 113°. The distribution of the magne- 
tization in this case is represented in the Fig. [l] (inset) . 
Extending this argument to the region M ~ 1, one finds 
the Doring result 7 = arccos(— 11/29) ~ 112°, showing 
that the rotation angle 7 must actually be a function 
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FIG. 2. (Color online) Bloch point for the maximum energy 
structure, 7 = 0. (top) Nondimensional magnetization inten- 
sity profile for v = 1, 10, and its angular distribution (inset), 
(bottom) Variation of the normalized slope at the origin as a 
function of v, and its asymptote Coo (dashed line). 



of the radial coordinate, although its slow variation con- 
firms the approximated validity of the ansatz ( 10 1 for this 
structure. 



Second, we consider case (ii), the global solution with 
sin 7 = 0, corresponding to a maximum of the energy 
structure. In this case we note that Eqs. (19 1, (20), and 



( 22 1 hold identically, and then we are left to solve 



DM(r) + {2v - l)M(r) - 2vM{rf = (29) 

where the last term comes from F'(r) = —M(r) in the 
magnetized domain r < R. The behavior of M(r) for 
7 = 0, 7r, shown in Fig. [2] is qualitatively similar to the 



M s {v) = y/(2v-l)/2v, 



(30) 



smaller than the v — > oo limit value of 1, relevant in 
the minimum energy case, where we neglected the de- 
magnetizing field. It is worth noting, that in spite of 
the similarity between ( 24 ) and ( 29 1 , the two cases are 



completely different: the global solution provide an exact 
solution of the magnetization profile in the spherical re- 
gion. It is not possible to connect the two solutions, the 
local solution corresponds to the minimum of 5(7) while 
global case to its maximum. Incidentally, we remark that 



the second (12 1 condition is violated for v < 1/2, show- 
ing that for small v the radial Bloch point does not exist 
as a stationary state. Near the origin the magnetization 
is linear with a slope C(u), represented in Fig. [2] that 
tends to C 00 {2v) 1 / 2 for large v. This slope determines 
the typical size of the singular core, 



i s = i/C{y) 



(31) 



which tends to zero when v — > 00, showing the singular 
behavior near the origin in the micromagnetic approxi- 
mation. The bottom panel of Fig. [2] shows the slope C, 
whose inverse det ermines the characteristic length of the 
singular core (31 1, normalized to (2v) x l 2 . 

In conclusion, we have shown that the interplay of 
exchange, dipolar and Landau energies, determines the 
internal structure of the Bloch point, beyond the mi- 
cromagnetic approximation. We have generalized the 
pure radial solution of Galkinap' and the constant mag- 
netization modulus solutions derived from the original 
Feldtkellei 6 ansatz. Noteworthy, the actual structure of 
a Bloch point eventually depends on the boundary con- 
ditions, that will determine the demagnetizing field and, 
near the singularity, select the rotation angle 7 (in gen- 
eral a function of the radius) An estimation of the size 
of the inner core gives £ s w 2 1 / 2 £ ~ 1 nm (for large 
v and typical permalloy parameters values, I sa 10^o), 
shows that in order to resolve the dynamics it is necessary 
to introduce quantum effects at least in a semiclassi- 
cal approximation. In fact, such effects will control the 
dissipation processes that lead to the actual topological 
transition. 
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